Diffusive-ballistic Crossover in ID Quantum Walks 
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We show that particle transport, as characterized by the equilibrium mean square displacement, in 
a uniform, quantum multi-baker map, is generically ballistic in the long time limit, for any fixed value 
of Planck's constant. However, for fixed times, the semi-classical limit leads to diffusion. Random 
matrix theory provides explicit analytical predictions for the mean square displacement of a particle 
in the system. These results exhibit a crossover from diffusive to ballistic motion, with crossover 
time on the order of the inverse of Planck's constant. This diffusive behavior is a property of the 
equilibrium average and does not require further interactions of the system with the environment. 
We expect that, for a large class of ID quantum random walks similar to the quantum multi-baker, 
a sufficient condition for diffusion in the semi-classical limit is classically chaotic dynamics in each 
cell. These results describe an interesting generalization of known quantum random walks and may 
have applications for quantum computation. 
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Recent results for non-equilibrium and transport prop- 
erties of extended classical systems with microscopic 
chaos P, 13 have suggested that one might explore the 
transport behavior of quantum versions of simple, ex- 
tended classical systems. One such system, studied here, 
is the multi-baker model for deterministic diffusion in 
one-dimension. In paper Q we introduced a quantum 
version of the classical multi-baker model as a combina- 
tion of local quantum baker dynamics with a transport 
of wave functions to neighboring cells. Rs simplest case 
is an example of the quantum random walks that have 
been considered previously (the Hadamard walk) 0, |^ , 
but other cases form a much wider class of ID quan- 
tum random walks. Quantum random walks have be- 
come of considerable interest over the past few years since 
they might be useful for implementing quantum random 
search algorithms if quantum computation becomes pos- 
sible. Quantum multi-baker models might be excellent 
candidates for implementation as part of current efforts 
to provide techniques and algorithms designed to take 
advantage of the possibilities inherent in quantum com- 
putation. Its experimental realization eventually may be 
feasible since both quantum baker maps 0, Q , and the 
Hadamard walks 5] , are now within experimental reach. 
The general case we consider allows a larger number of 
quantum states participating in the random walk, and 
widens the range of possible physical applications. 

In this Letter we report on the behavior of the equi- 
librium mean square displacement (m.s.d.) for a quan- 
tum particle whose classical dynamics is governed by 
the deterministic, diffusive multi-baker process. Since 
one-dimensional processes in extended quantum systems 
show ballistic motion, for translationally invariant sys- 



tems, or localization, for disordered ones, it is interesting 
to study how each system makes the transition to diffu- 
sive behavior in the classical limit. Here we concentrate 
on the translationally invariant case. A natural question 
is whether or not an interaction with the external world, 
or decoherence, is necessary to eventually restore diffu- 
sive motion in the classical limit. Our work here shows 
that this is not the case for the m.s.d.: for systems that 
we study there is a crossover from short-time diffusive 
behavior to ballistic motion in the long time limit. The 
crossover time is on the order of the inverse Planck con- 
stant, i.e. the Heisenberg time, rather than the Ehren- 
fest time, on the order of the log of Planck's constant. 
The chaotic classical motion of our system allows the 
use of random matrix theory (RMT), which leads to ex- 
plicit expressions for the mean square displacement of a 
quantum particle. The comparison of the RMT analytic 
results with results of numerical evaluation of the exact 
formula for particular systems is very good with some 
interesting exceptions. Our calculation has a classical 
counterpart 2] which has the same result for the m.s.d. 
as obtained here in the semi-classical limit. 

We begin with the classical multi-baker map [1, ■ It 
is a deterministic model of the random walk on a one- 
dimensional lattice, where the direction of the jump is 
determined by the internal state of the particle. We use 
the baker map as a model of internal dynamics. That is, 
the multi-baker map, M, is a composition of two maps 
M = B oT. First, the phase points are transported to 
neighboring cells according to T(n,x,y) = (n + l,x,y), 
for < a; < 1/2, and (n-1, a;, y) for 1/2 < x < 1. Then a 
baker map B acts on the x, y coordinates of each cell, n, 
separately, according to B{n, x, y) = (n, 2x, J//2), for < 



2 



X < 1/2, and (n,2a;-l,(l + y)/2), for 1/2 < x < 1. The 
combination of these two maps is the multi-baker map 
which is a time-reversible, measure preserving, chaotic 
transformation, with evolution law 



M{n,x,y) 



(n + l,2a;,f). 



for < a; < 1/2, 



(n-l,2a;-l,i±^), forl/2<a;<l. 



It is the simplest area-preserving model of a random walk. 

To quantize the dynamics we first quantize i? in a 
single cell. The horizontal direction of the torus [0, 1)^ 
is taken to be the coordinate axis, while vertical axis 
corresponds to the momentum direction. To obtain 
the Hilbert space |0, 0, we take the subspace of 
the wave functions on a line whose probability densi- 
ties, are periodic in both position and 
momentum representations, respectively: ^'(a:: -I- 1) = 
exp{i2TTipq)^{x), ^{p -I- 1) = exp(z27r(/3p)^'(p), where 
ipq,ipp G [0,1) are phases parameterizing quantization. 
The quantization of the baker map requires the phase 
space volume to be an integer multiple of the quantum of 
action pH Ull , so the effective Planck constant must 
be h = 1/N, where N is the dimension of the Hilbert 
space. The space and momentum representations are 
connected by a discrete Fourier transform G^i'-Pq, fp)-— 
{Pk\qi} — N^^/^ exp{~i2TTNpkqi)- The discrete positions 
and momenta are qi = {I + ipq)/N, pk = {k + ipp)/N. 
The Hilbert space H can be decomposed into "left" and 
"right" subspace H = Hl ® Hr, and "bottom" / "top" 
spaces H = Hb ® Ht, which are N/2 dimensional, 

(^Hl when {qiYi) = for ; = N/2, . . . , N -I, e Hb 
when (pfcl^-) = for fc = N/2, . . . , TV - 1. 

Having constructed the Hilbert space one looks for a 
family of unitary propagators parameterized by = 1/ft. 
which go over into the classical map in semi-classical 
limit. Technically, one requires the Egorov condition to 
be satisfied, which means semi-classical commutation of 
the quantum and classical evolution 1121. The unitary 
operator for the quantum baker map |lllll2l is given 



by Un := G 



N 



for even TV. Other exam- 



Gjv/2 

Gn/2 

pies and discussions of issues concerning the quantization 
of area-preserving maps can be found e.g. in (lij . 

Since the quantum multi-baker is a model of particle 
with N internal states jumping over lattice of length L, 
we use product Hilbert space H = ® to describe 
it. The dynamics is implemented in two steps. First one 
shifts states from right subspace at cell n {Hji{n))to right 
states at cell n — 1, and states HL{n) into left states at 
cell n + 1, which gives the quantum transport operator 
T . Then on each of the cells one acts independently with 
a quantum baker operator Um, which corresponds to the 
classical map B. Let us write the states \^) G H in the 
position basis X]n=o Sj=o^ j)- Then the only 

non-zero matrix elements of the quantum multi-baker op- 
erator M are of the form {n + 1, j\M\n, k) = {j\UN\k), 
for k e -ffi,(n), or {n - l,j\M\n,k) = {j\UN\k), for 



k G H[i{n). For the translationally invariant case studied 
here, the m.s.d. can be expressed entirely in terms of the 
properties of the baker operator, B, cf. Eq. (2), and we 
do not need to make explicit use of the structure of M. 

The central quantity of interest, is the expression for 
the mean square displacement (m.s.d.) of a quantum 
particle in the chain, defined as the average value of 
the mean square displacement taken with an equilibrium 
density matrix for the system. This is a uniform distri- 
bution of probabilities along the chain, geq = Inl/NL; 
{A) := Tr {geq A) = Tr {A)/NL, L is the length of 
the chain (we assume periodic boundary conditions and 
where N is the dimension of the Hilbert space. Then 
the m.s.d. is simply {{Arfit)) = {{M^^^rM* - rf). If 
we define the velocity operator v := M'^rM — r then the 
m.s.d. can be written as ((Ar)2(t)) = X]m,n=o(^m"")' 
where u„ := M^"wM". Time invariance of gieq im- 
plies (vmVn) — {v„i-nVo)- Thus wc Can express the 
m.s.d. in terms of the velocity autocorrelation function 
Cn = {v„vo): 



t-i 



{{Arf{t))^t{v')+2j2(.t-n)C„ 



(1) 



We use a "coarse" position operator r defined by 
r\n, k) := n\n, k). We calculate the coarse velocity oper- 
ator on the line using the definition given above, obtain- 
ing [ill, v{n,k;m,l) = ±Sk,iSn,m, with + for / < N/2, 
— for other /, and then put it on the circle to enforce 
translational invariance. This form can be understood by 
observing that for the translationally invariant case, the 
coarse velocity ±1 denotes a translation of the quantum 
state one cell to the right or left. An identical form for 
the coarse velocity occurs in the classical multi-baker as 
well 0. Then the velocity autocorrelation function can 
be reduced to the trace over states in a single cell (1^] 

Cn = — Tr [M^"i>M"v] = — Tr [Bt"JB"J], (2) 
LN N 

where we note that J is the velocity operator v reduced 
to a single cell. In the position representation J is given 



Assuming the properties of the 



local propagator, B, are known, one can express the ve- 
locity autocorrelation function in terms of its spectrum 
and eigenvectors, where B\k) = cxp(z<y9fc)|A;). It immedi- 
ately follows that 



Gn — 



1 



(3) 



where Jjk '■— {j\J\k). Since J has a very simple form, 
one sees that 



\J,k\' = Tr J^ = N. 



(4) 
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Substituting these results in formula we obtain 
1 .'^r-, , 1 



33 I 



TV 



(5) 



3^k 

i 



n=l 

-E 



2 „;„2 {ajkt) 



<^7"fc 



TV 



,2 ajfc 



(6) 



where Ojk '■= fj — Vfc- Whenever there is a degeneracy, 
ctjk — 0, we replace the ratio of sines by t^. In certain 
cases it is possible that the coefficient of the ballistic, 
i^, term may be zero; see Figure |21b and the following 
discussion. 

We now evaluate the m.s.d., Eq. jSjl using random ma- 
trix theory and compare the results with numerical eval- 
uations. To apply RMT we consider the velocity auto- 
correlation function C„ given by Eq. |(2Jl, and separate 
the terms on the right hand side into diagonal, j — k, 
and non-diagonal, j ^ k, terms. We suppose that B is 
drawn randomly from either the COE or CUE ensembles, 
although the numerical results present a more general be- 
havior. We assume the distribution of matrix elements is 
independent of the distribution of elements of eigenvec- 
tors (see section 8.2 of Il6|, and Using TrJ^ = N, 
one sees that the ensemble average of the mean square 
displacement takes the form 



(((Ar)2)> 



t + t{t 
+2{N- 



t-i 

l)U^k\')Y.{t-n){e^n- (7) 



We replace the matrix elements \Jjj\'^ and |Jjfep by 
their average values and (|JjfcP), respectively. 

Straightforward calculation gives = k/{N + 

k), where fc = 1 for CUE, and 2 for COE. Averaging 
Eq. Q, we obtain (|JjjP) + (N - l)(\Jj^k\'^) = 1 and 
thus {\Jj^k\^) = N/[{N + k){N - 1)]. Then we need 
to calculate the average value of the exponential factor 
exp[i((/jj —ipk)n]. For this calculation we need the expres- 
sion for the pair correlation function R{^j, (fk) in the two 
ensembles, so that we can express the average value as 

r2TT p2tt 



[i(yi-Vfc)n]^ _ 



N{N-1)' 



These correlation functions are given in the litera- 
ture For the CUE, one finds that 



47r2 



sm 



2 N{ipj-ipk} 



N'^ sin 



2 iVj-Vk) 



A straightforward calculation leads to 



1 



,[ii'Pj-Vk)n]^ ^ 







for n = 
^ for n < N, 
for n > TV. 



(8) 



Using this estimate for the average the result, we find 
that the m.s.d. is given by 



((Ar)^(i)) 



t + 

k 



N+k 



«(*-!) 
N+k 



t^ + f - 



1 

N{k_ 



t-2 

3(7V-1) 

'1) 



3(Ar-l-fc) 



for t < N, 
for t > N. 



Note that the "super-ballistic" term only occurs for 
t < N, where it is typically less than or on the order of 
the linear term, t. The exact result for the COE ensem- 
ble has a correction arising from an additional term in 
the pair correlation function. This correction is rather 
lengthy to write and is negligible for both short and very 
long times, with the maximum deviation of at most five 
percent occurring at t = TV. The details will be given 
elsewhere Figure ^ shows the estimates for the two 
ensembles for TV ~ 200. 
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FIG. 1: Log-log plot of the ensemble averages of the mean 
square displacement using RMT. The COE results are the 
two close curves, where the lower is the result given in Eq. © 
for k = 2, while the higher curve was obtained using the 
full pair correlation function for COE. Three asymptotic es- 
timates t,t'^ /N,2t^ /N are also plotted. The inset shows the 
region t = 100 to t = 300 where the differences between the 
two COE results are most pronounced. 



We compare these predictions with numerical results. 
For almost every choice of phases ipq,fp defining the 
quantization, the evaluated formula ((HJ gives results be- 
tween the COE and CUE average predictions for TV 
greater than 100. The Balazs-Voros phases {(fq = ipp = 
0) [Isj yield exceptionally good agreement with CUE av- 
erage for all values of TV. Figure (|21a) compares the exact 
result with the RMT averages in this case. 

An interesting exception to these results occurs when 
fn-\- — 1, e.g. for the Saraceno case {^pq — ipp = 1/2) 
[ll|. For these special values of the phases, the proba- 
bility of finding the system on each half is 0.5 for every 
eigenstate, so that Jjj = 0, and there are no degeneracies. 
Thus for long times the m.s.d. oscillates around the time 
average value, ^j^k |^jfcP/(2TVsin^ after an initial 
diffusive growth. To see how such oscillations are possi- 
ble, consider a more general class of systems where the 
local quantum baker map is replaced by an arbitrary uni- 
tary operator B ^| . For this larger class of systems one 
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FIG. 2: Comparison of RMT estimates with the numerical 
evaluation of the formula Q for the m.s.d. in case of quantum 
multi-bakers with a) Balazs-Voros phases (generic case), b) 
Saraceno phases (exceptional, localized case). Plots are in 
double logarithmic scale, the inset in figure (b) shows the 
results for much longer time in normal scale. 



easily proves that < ((Ar)^) < i^, both classically and 
quantum-mechanicaUy, and it is easy to construct exam- 
ples realizing both extrema. Moreover, if i? is a right-left 
exchange operator, which merely replaces left and right 
states, the m.s.d. oscillates between and 1. Similar phe- 
nomena over much longer time-scales may take place in 
the case discussed above. Figure (|21b) shows reasonable 
agreement with the short-time classical diffusive behav- 
ior also in the case of Saraceno quantization. We expect 
that for large N the phases should not matter for times 
up to the Heisenberg time. 

Summary: We have considered the mean square dis- 
placement of a particle whose dynamics is governed by 
the propagator for a quantum multi-baker map. After 
deriving a general formula for the m.s.d., we evaluated it 
by means of RMT and also by numerical methods. Ran- 
dom matrix theory provides explicit, analytic expressions 
for the m.s.d., as well as for the velocity autocorrelation 
functions that determine it. Our study shows that these 
analytic expressions depend somewhat on the nature of 
the circular ensembles used, either unitary or orthogonal. 
Comparison with numerical results shows that neither 
of the ensembles is superior: depending on quantization 
phases one or the other gives a better representation of 
the data, but usually the experimental curve lies between 
the two predictions with numerical ballistic coefficient 
having values between and for generic (bal- 



listic) systems. The analytic expressions for the m.s.d. 
allow us to study the transition to classical behavior in 
detail. We find that, on the average, there is a smooth 
transition from quantum to classical behavior as Planck's 
constant approaches zero, and that for non-zero values of 
h the classical behavior persists up to times on the or- 
der of , and we see no need for further interactions 
with the environment, for a well behaved classical limit 
for the models we study. Thus the analytic expressions 
for the velocity autocorrelation functions and m.s.d. pro- 
vide a powerful tool for studying the quantum-classical 
transition for simple extended systems. Given that it is 
possible to realize a quantum baker map experimentally, 
the quantum multi-baker map eventually may have ap- 
plications for quantum computations. 

We show elsewhere that most of the results presented 
here are valid for a much wider class of quantum sys- 
tems 0. 
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